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Abstract 

We develop a semi-classical approximation for the scar function in the Weyl-Wigner represen- 
tation in the neighborhood of a classically unstable periodic orbit of chaotic two dimensional 
systems. The prediction of hyperbolic fringes, asymptotic to the stable and unstable manifolds, 
is verified computationally for a (linear) cat map, after the theory is adapted to a discrete phase 
space appropriate to a quantized torus. Characteristic fringe patterns can be distinguished even for 
quasi-energies where the fixed point is not Bohr-quantized. Also the patterns are highly localized 
in the neighborhood of the periodic orbit and along its stable and unstable manifolds without any 
long distance patterns that appears for the case of the spectral Wigner function. 
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I. INTRODUCTION 



The Gutzwiller trace formula provides a tool for the semi-classical evaluation of the energy 
spectrum in terms of the periodic orbits of the system. However the number of long periodic 
orbits required to resolve the spectrum increases exponentially with the Heisenberg time T H 

Q. 

The semi-classical theory of short periodic orbits developed by Vergini and co-workers 
[ll is a formalism that allows to obtain all the quantum information of a chaotic Hamiltonian 
system in terms of a very small number of short periodic orbits. In this context, the scar 
functions play a crucial role. These wavef unctions, that have been introduced in several 
previous works [3]- [3], live in the neighborhood of the classical trajectories, resembling the 
hyperbolic structure of the phase space in their immediate vicinity. This property makes 
them extremely suitable for investigating chaotic eigenfunctions. Recently it has been shown 
that the matrix elements between scar functions provide information about the heteroclinic 
classical structure {(], Q] • 

In order to perform further developments to this semi-classical theory of short periodic 
orbits, it is important to provide a general expression in phase space for the scar functions 
explicitly in terms of the classical invariants that generates the dynamics of the system. This 
is the purpose of this paper. In this context the choice of the Weyl-Wigner representation 
provides a phase space vision which best allows a quantum classical comparison [13I ]. 

The semiclassical expression here deduced predicts characteristic hyperbolic patterns lo- 
cated in the neighborhood of the periodic orbits and along its stable and unstable manifolds. 

In spite to compare the general expression here founded with a "realistic" system the 
cat map is employed, i.e. the quantization of linear symplectic maps on the torus. As was 
shown by Keating 21] in this case the semiclassical theory is exact, making this maps an 
ideal probe for our expressions. After the formalism is adapted for a torus phase space, we see 
an important agreement between the semi-classical construction of the and the numerically 
computed Wigner scar functions for the cat maps. The characteristic hyperbolic patterns 
are clearly discernible even for non Bohr-quantized values of the classical action. 

Although most studies in phase space have adopted the Hussimi representation 
this must be interpreted as a Gaussian Smoothing of the Wigner function in an area of size 
h that usually dampens the fine structure extending along the stable and unstable manifolds 
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This work is a first step to further objectives that includes the investigation of the matrix 
elements between scar functions. 

Section 2 deals with the definition of scar functions and develops its relationship with the 
spectral operator. Then, the Weyl-Wigner representation of the scar function, here called 
Scar Wigner Function, is studied. Mainly, its local hyperbolic form of the neighborhood of 
a single periodic orbit is obtained only in terms of classical objects. Also, this expression is 
compared to the already studied spectral Wigner function [fjj] . 

Section 3 is devoted to study the particular case of the cat map where not only the semi- 
classical theory is exact but also the linear approximation is valid throughout the torus. The 
scar Wigner function are then clearly visualized as a hyperbolic fringe pattern that agree 
with the semiclassical expression derived in section 2. It has to be noted that we here choose 
values of N such that the scar states are not eigenfunctions of the cat map as Faure et al. 
[12I ] shows it happens. Also the fringe patterns here analyzed lie entirely within the scope 
of standard semiclassical theory, restricted to fairly short times. The interesting question 
concerning homoclinic recurrences depends on dynamics for longer times. 



II. SCAR WIGNER FUNCTIONS 



Scar function states are the main object of study of the current work. According to 
the scar function \<px,<t>) of parameter constructed on a single periodic point X = 



(P, Q) is defined as 



T 

<Px*) = J^dte* cos (fyv*\X) (1) 



where T = lnh is the Ehrenfest time, \X) is a coherent state centered in the point X on 
the periodic orbit and IP is the unitary propagator that governs the quantum evolution of 
the system. This wavefunctions have been shown, in the Hussimi representation, to live in 
the neighborhood of the trajectory, resembling the hyperbolic structure of the phase space 
in their immediate vicinity [121 ]. 

On the other hand the spectral operator (or energy Green function) is defined as 



00 
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G Et£ = / dte^lPe-^ (2) 
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with U — e a m and 



j- so that 



Ge,e — — r 
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H-E) +e 



M-E-ie, 

5 e (E-H^)=J2 IV>n> (Vnl ^ - 



(3) 
(4) 



where |^ n ) are the Hamiltonian eigenfunctions with energy E n and 5 £ (x) is a normalized 
function whose width e can be taken to be arbitrary small, so that 5 e tend to the Dirac 5 
function as e — > 0. Then Ge, £ is not a pure state but a statistical mixture. In the limit e — » 
the spectral operator is a comb of delta functions on the eigenangles, whose amplitudes are 
the corresponding individual density matrices \i/; n >< ij) n \. For values of e larger than the 
mean level spacing, several eigenstates contribute to the spectral operator in a Lorentzian 
-like smoothing of energy width e. 

It is important to mention here that the cosine term in ([1]) has been chosen so that the 
energy spreading is minimum 7| but can be replaced by any damping term with charac- 
teristic time T as for example, e~~^ with e ~ h/T. So that, it is possible to establish a 
relationship of the scar functions with the spectral operator 



\fx,4> >— Gtyfi/TiX > 



(5) 



that is, the scar function is the spectral operator acting on the coherent state centered at 
the point X on the periodic orbit. 

The purpose of this work, is to study the scar function in phase space by means of the 
Weyl Wigner representation p x ,(x), here call the Scar Wigner function, 



PxA x ) = tr R xPx, 



(6) 



where px,^ = \(px,<t) (<Px,<i>\, is the density matrix of the scar function and R x are the set 
of reflection operators thought points x = (p,q) in phase space [l5, Q. The Weyl- Wigner 
representation, defined through the set of reflection operators (see Appendix), has the ad- 



vantage to show structures of size lower than h 



a Gaussian smoothing on a region of size h 
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171 ] while the the Hussimi representation is 
and hence has a lower resolution. 



For the case of the spectral operator its Wigner function 



G E ,e(x) = W(E,e,x) = tr 



(7) 
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is known as the Spectral Wigner function, first introduced by Berry [19] , has been recently 
shown to present important scarring features pjj some of which would be present in p x , (x) 
as will be seen in this work. 

It must be noted that by construction 



px,<p = \<px,4> >< ( Px,4 



is a pure state while Ge, £ is a statistical average but both are related by 



(8) 



Px,<t> — G h ,h/t\X >< X\G^ h/T , 



then for the scar Wigner function 



PxA X ) = tr R xPX,, 



< X \ G ^,n/T^ G ruj>,n/T\ x > 



(9) 



Inserting the definition of scar function ([I]) and (|HJ) in it can be seen that 



pxA x ) 



T T 
2 I 2 



dt'dte m cos I -- I e~ iq>v cos I '— I < X\U~ V BJP\X > 



T j T_ 
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Tit 



T 



-i<j>t' 



irt' 



T 



t'n frt\ 



(10) 



Let us now use the decomposition of the propagator in terms of reflection operators 15] 

L 



-U / dxU\x)R x 



(11) 



where j dx is an integral over the whole phase space and L is the number of degrees of 
freedom. So that, 



(12) 



The coherent states on points X = (P, Q) in phase space are obtained by translating to X 
the ground state of the harmonic oscillator, its position representation is 

"mui\ i 



(q W = ( 



exp 



2h iq - Q) +l h[ q -2 



(13) 



For simplicity, unit frequency (u — 1) and mass (m = 1) are chosen for the harmonic oscil 
lator without loss of generality. The overlap of two coherent states is then 

{X-X'f 



{X \X') = exp 



Ah 



2h 



X AX' 



(14) 
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Where the wedge product 



A A X' = PQ' - QP' = (JX) .X', 

with the second equation also defines the symplectic matrix J . As is shown in the Appendix 
the action of the reflection operator R x on a coherent state \X) is the x reflected coherent 
state 

R x \X) = exp {^-X A x^j \2x - X) (15) 
and the product of three reflections also gives a reflection 



R X2 R X R X1 = e i A3 ^ Xl ^R 



X2-X+XX 



(16) 



where A 3 (x2, x±, x) is the area of the oriented triangle whose sides are centered on the points 
x 2 ,Xi and x respectively. 

Then applying (TIB"]) then (fT5|) and performing the overlap (fl^j) it is shown that 

,2~ 



(X\R X2 R X R X1 | A) = e* A3 ^^^exp. 



(X-x R ) 

h 



with 



x R = x 2 - x + x 1 . 

Inserting (TTT|) in f[T^]) in order to perform the double phase space integrals 



(17) 



(18) 



2L 



1 = [ ^ ) / dx ± / dx 2 U~ t '(x 2 )U t (x 1 ) (X\ R X2 R X R X1 \X) 



1 



2L 



dx x / dx 2 U- t ' '{x^U'ix^e^ 2 ^ 



exp 



(a - x R y 
h 



(19) 



If the exponential term is omitted in ( 1191 . the double phase space integral is simply performed 
using the product of symbols in phase space [l5] that is 

(x) = U-*U*(x) = (J^j J dx x J dx 2 U-* (x 2 )f/ t (x 1 )e^ A3(x2 ' a;i ' a;) (20) 

but the presence of the exponential term has to be kept into account so that 

(A - x R {x)f 



I = u*~* ( x ) exp 



h 



(21) 
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where now Xr(x) is a point in phase space that only depends only on x. Another approach 
to obtain this result is to perform in ffT9l) the semi-classical approximation for the propagator 



U'.Jx) 



E 



det M* + l 



-exp l -Sjt{ r 



(22) 



where the sum is performed over all the classical orbits 7 whose center lies on the point 
x 15]. Then S^t(x) is the classical center generating function of the orbit, M* = 9 



stand for the monodromy matrix and a* its Maslov index. 
Inserting (1221 in f|T9|) . we get 



EE< 

71 72 



cfxi / fix 



exp | [S^(xi) - S^(x 2 ) + A 3 (x 2 ,Xi, 



X 



|det M"/ 



1| 5 [detM^ + 1| 5 



exp 



- x R y 
h 

(23) 



Now the integrals are performed using the stationary phase approximation. In this case the 
phase <fi(x 2 ,x 1 ,x) = S^(xi) - S^(x 2 ) + A 3 (x 2 ,Xi,x) is stationary (i.e.) 

<90(x 2 , Xi, x) 



dxi 
d(j)(x 2 , Xi, x) 
dxo 








(24) 
(25) 



for x 2 (x), x\(x) such that the canonical transformation generated by (xi) is combined with 
that generated by Si ,(^2) to result in a new canonical transformation generated by S7_ t ,(x). 



That is, the canonical transformation corresponding to the trajectory 7* * is obtain as the 
co.npos.tion of the canon.ca. transformation ofthe orbits 7 < and respectively. For c.arity 
this situation is depicted in Fig 1. (see also [15l|). 
Hence 



E- 



ex P i[S*-*'(x)] 



IdetAft-f 



1 



exp 



(X - x R (x)Y 

h 



(26) 



with Xr(x) = x 2 (x) — x + Xi(x) corresponding to the end point of the orbit 7* or to the 
initial point of 7^* as shown in fig. 1. Denoting x_ the initial point of the orbit 7J, its time 
t evolution is xr(x) 

x.(t)=x R (x). (27) 
While I. is the time t — if evolution of x_ 



X-(t — t') = x + , 



7 




FIG. 1: In order to correctly compose the canonical transformations the end point of the orbit 7*, centered 
in Xi and the initial point of 7^* , centered in x^ must coincide in xr. The resulting classical orbit 7*"* 
that joins the points X- and x+ and has center point x is the composition of the orbits 7* and 73"* . 

and x is center point of x_ and x + . That is, 

x = X.+X+ = x^+x_(t-t') ^ 

The resulting integral (126]) defined in ffl9l) is then inserted in (fT2l to obtain the semiclas- 
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sical approximation of the scar Wigner function 



T T 

py^x) = lljl dtdteWt-e) cos ^ cos (^L 



„sc 



x > e 7 ^- L - z jexp 

; IdetM*-*' + 



2 



(X-x fl (o;)) 2 
ft 



(29) 
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So as to study p^f^x) i n the neighborhood of the periodic point X of the map, the main 
contribution in the sum over classical orbits in (I2"9"j) will come from the particular periodic 
orbit that lies the closest from the periodic point X. Other orbits contribution will be highly 
damped by the exponential term involving X — xr(x). Then, only this particular orbit will 
be taken into account. 

For this purpose let us define x' = x — X. In the same way, x_ the initial point of the 
classical orbit 7*"* can be written as x_ = X + <5_ so that its time t evolution 

x R = x-{t) = A^x_ = (X + <L) = X + M*<L. (30) 

Where in the last expressions the flux in the neighborhood of X has been linearized and 
M* is the symplectic matrix denoting this linearized time evolution. Note that for the fixed 
point M*(X) = X. Also for the center point defined in ff28l) 

X = X + S= x - +x f-V = X + l - +X 2 + M '<" S - = X + (M^ + l) S f. (3D 
Inverting this last expression 

5- = 2 [M 1 - 1 ' + \y l x' (32) 

that is inserted in (!30~!) to obtain 

^= x+2 (ssrTi) x ' (33) 

and finally, 

M t 

X - x R = -2 , - J rV. (34) 

The eigenvalues of the symplectic matrix M* are exp(— At) and exp(At), (A is the stability 
exponent of the orbit) corresponding to the stable and unstable directions generated by 
the vectors £ s and £ u respectively. Let us define q' and p' as canonical coordinates along 
the stable and unstable directions respectively such that x' = (p 1 , q') = q'£ s + with 
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( M A( S = 1. Then, using fl34j) with the diagonal representation of the symplectic matrix, the 
scalar product (X — xr) 2 = (X — xr) . (X — xr) takes the form 

1 



(X - X R ) 



cosh 2 (A*^) 

where £ 2 = £ u .£ u and ^ = £ s .£ s . It can also be easily shown that 



(35) 



| del M* + 1 1 f = 2 cosh j ^ 



The center action for the orbit that is close to the periodic orbit is given by 



Sl(x) =tS x +x'B*x' + 0(x' 3 ) 



(36) 



(37) 



where Sx is the action of the periodic orbit (fixed point) for which the Maslov index at = ia 7 . 
Let us define the action Sx = Sx + fa* 7 in order to include the Maslov index in the action. 
B l is the symmetric matrix such that 



. I- Mi 



' 1 - M* 



with 



J 



-1 

1 



(38) 



(39) 



Thus B l is the Cayley parametrization of Ml. Using the stable and unstable directions as 



coordinate axes, the result is that 



JBt 



tanh(iA/2) 

-tanh(tA/2) 



(40) 



Inserting (jJUJ) in fl37j) and putting together with fl36l) and fl35|) in ff29l we obtain that the 
contribution to the semi-classical scar Wigner function from the orbit that lies close to the 
periodic orbit X takes the form 

T T 



dt'dt exp % { I <f) + ^ ) (t - t') - ~p'q' tanh 



(t - f) x 



cos (f ) COS (f^) 



2 cosh 



x exp 



(t-t')A 
2 



(41) 



cosh 2 (A^) 

this last expression shows that the dependence of the scar Wigner function on the phase 
space variables has two aspects, a phase oscillating term only depending on the product 
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p'q' and a damping term. The phase oscillating terms shows phase coherence along the 
stable and unstable directions where the damping factor is null. Phase coherence also 
holds where the product p'q' is constant, i.e. along each successive hyperbola that has the 
stable and unstable manifolds as asymptotes. This implies in phase oscillations across the 
hyperbolae, the amplitud of the oscillations decrease with increasing A and will be maximal 
for (f) corresponding to the Bohr-quantized orbit. While the damping term, implies that away 
from the asymptotes the amplitude of the phase oscillations presents a Gaussian decreasing 
with increasing the distance. 

This amplitude decaying is not present in the Spectral Wigner function that has strong 
oscillations away from the periodic orbit and its stable and unstable manifolds [ljj]. This 
facts implies in important superpositions of the hyperbolic patterns from different orbits 
that usually washes out all the structure. 



III. SCAR WIGNER FUNCTIONS FOR THE CAT MAP 



Now the present theory is applied to the cat map i.e. the linear automorphism on the 
2-torus generated by the 2x2 symplectic matrix Ai, that takes a point X- to a point x + : 
x + = M.x^ mod(l). In other words, there exists an integer 2-dimensional vector m such 
that x + = M.X- — m. Equivalently, the map can also be studied in terms of the center 
generating function [22]. This is defined in terms of center points 

X + + £_ 



X 



(42) 



and chords 



where 



X j 



x_ 



-J 



dS(x, m) 
dx 



(43) 



S(x, m) = xBx + x(B — J)m + -m(B + J)m 



(44) 



is the center generating function. Here B is a symmetric matrix (the Cayley parameterization 
of M, as in (gO])), while 

n 1 

(45) 








1 


J = 








1 
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We will study here the cat map with the symplectic matrix 



M 



2 3 
1 2 



and symmetric matrix B 



1 



(46) 



This map is known to be chaotic, (ergodic and mixing) as all its periodic orbits are hyper- 
bolic. The periodic points xi of integer period / are labeled by the winding numbers m, so 
that 

Pi 
Qi 



Xi 



r 



m. 



(47) 



The first periodic points of the map are the fixed points at (0, 0) and (|, |) and the periodic 

■> i ■ i 2" '!' L l 2'(l' ' \2> £i)]' e) ' §)]' [(^' 5)' (2' an d 



orbits of period 2 are [(0, §) , (|,0)], [(§, \) 
[(0, |), (0, §)]. The eigenvalues of M are e" A and e A with A = ln(2 + \/3) « 1.317. This 
is then the stability exponent for the fixed points, whereas the exponents must be doubled 
for orbits of period 2. All the eigenvectors have directions £ s = (— |) and £ u = (1, 
corresponding to the stable and unstable directions respectively. 

Quantum mechanics on the torus, implies a finite Hilbert space of dimension iV = 
and that positions and momenta are defined to have discrete values in a lattice of separation 



jj [l6j , |2(| . The cat map was originally quantized by Hannay and Berry 20| in the coordinate 
representation: 



< q* U M a< >= — 



N V 



jk + f 



(48) 



where the states < g|qj > are periodic combs of Dirac delta distributions at positions 
q = j/Nmod(l), with j integer in [0, N — 1]. In the Weyl representation 



map has been obtained in 



16| . the quantum 



221 as 



Vm{x) = i y e *«PW 

|det(A4 + l)|5 V 
|det(A< + 

where the center points are represented by x = (77, ^) where a and b are integers in 
[0, N — l] for odd values of N [if]]. There exists an alternative definition of the torus Wigner 
function which also holds for even N. However, it is constructed on the quarter torus and 
this compactification scrambles the hyperbolic patterns. 
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The fact that the A4 matrix has equal diagonal elements implies in the time reversal 
symmetry and then the B matrix has no off-diagonal elements. This property will be valid 
for all the powers of the map and, using we can see that it implies in the quantum 
symmetry 

V l M (p,q) = {V l M (-P,Q)y = (UW, -<?))*■ (50) 

It has been shown [20] that the unitary propagator is periodic (nilpotent) in the sense 
that, for any value of N there is an integer k(N) such that U^ V ' ) = e 1 ^ . Hence the eigenvalues 
of the map lie on the k(N) possible sites 

i(2mTT + </>) 



exp 



k(N) 



1 < m < k{N). 



(51) 



For k(N) < N there are degeneracies and the spectrum does not behave as expected for 
chaotic quantum systems. In spite of the peculiarities in this spectra of quantum cat maps, it 
is likely that non-degenerate states are typical of chaotic maps, such as very weak nonlinear 
)erturbations of cat maps that are known to have non degenerate spectra |23J]. Eckhardt 



24| has argued that typically the eigenf unctions of cat maps are random. 



The Scar Wigner Function on the torus depends on the definition of the periodic coherent 
state 
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26l | . with < p >= P and < q >= Q. In accordance to ([13 



< X|q fc >= ^2 ex P \ 2 nN 



u + Q J IN? - ipu + 1 - «m 



j=-oo 

The Scar function is then defined on the torus as 

T/2 



(52) 



-T/2 ^ ' 



(53) 



Note that for maps, time only takes discrete values, then the time integrals in and 
1]) are in this case replaced by summations. Now, the scar Wigner function on the torus is 



pxA x ) 



Tr 



(54) 



where the trace is now taken on torus Hilbert space and R x are the periodic reflection 
operators on the torus 16]. 

In order to construct the semiclassical Scar Wigner Functions on the torus we have to 



periodized the construction. This is done merely using the recipe 



161 ] that for any operator 
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its Weyl representation on the torus A(x) is obtained from is analogue in the plane A(x) by 

OO OO / 7 ■ \ 

A(x)=J2 E (-l) 2ia+2kb+lkN A(x+ [ -^), (55) 

j=—oo k=— oo 

This leads to the property 

A(x + ^p-) = (-l)^ a+2kb +^ kN A(x), (56) 

for the torus Weyl symbol. In the case of the Scar Wigner function, the phase factor in 
fl56|) . leads to four images of the scar pattern, supported by the integer lattice. The images 
centered on (P, Q), (P + 1, Q) and (P, Q + \) are all identical, whereas (P + §, Q + 1) centers 
a pattern which is the negative of the other ones. This fact has already been studied for the 

r~i 

Wigner function of coherent states [14j . 

In figure 2 we compare the exact Scar Wigner function for a cat map with N = 223, a 
value for which the quantum map has no degeneracies. In figure 2 (a) with the semiclassical 
approximation, correspondingly for h = 1/ (2ttN), in (d), for the periodic point at (1/2, 1/2) 
whose action is Sx = 0.75 for a value of <fi that does not Bohr-quantize the orbit. Figure 
2.(b) and 2.(c) show respectively horizontal and vertical sections of the objects plotted in 
figures 2. (a) and 2.(d) in the neighborhood of the periodic point. 

It can easily be observed that the four images of the scar patters are present both for the 
exact and semiclassical scar Wigner functions. Also, the semiclassical approximation prop- 
erly describes the overall behavior of the exact dynamical system with detail of oscillations 
of the order 0.02. 

Although, deviations are present, particularly seen in figure 2.(b), they are due to the con- 
tribution to ( HIT) of longer orbits. In particular, the period-2 orbit on the points [(0, |)(0, |)] 
has (1/2, 1/2) as center point. This periodic orbits only have contributions for the adequate 
values of t — t', in the specified case t — t' must be odd for the period-2 orbit to contribute. 

IV. DISCUSSION 

As was already observed for the case of the spectral Wigner function, the imprint of the 
classical hyperbolicity on the Scar Wigner function is so clear that it can be detected even 
for quasi-energies that do not correspond to a Bohr-quantized periodic orbit. 
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The general features exhibited by our calculations should also be discernible for nonlinear 
systems as, indeed, our deduction was not restricted to cat maps. Though the theory in 
section II is only local, we conjecture that distorted hyperbolae asymptotic to curved stable 
and unstable manifolds will bear fringes reaching out from the periodic point. In the case 
of a chaotic Hamiltonian for two degrees of freedom, this pattern should emerge in two- 
dimensional sections cutting the periodic orbit at a point. This plane should be transverse 
to that of the orbit, and in the energy shell. 

The cut-off with time in our definition of the Scar Wigner function affords equal treatment 
to all periodic orbits in the denominator of ( 14T1) . The reason why the main contribution 
comes from only one hyperbolic orbit is that for an orbit of period n to have coherent 
contribution in (f¥T|) t — t' must have specific values that differ in n. That is, there are only 
— terms instead of 2T in (jHJ). Also, T increases only logarithmically with h so that the 
semiclassical limit A — > nX implies in a a cut off for longer orbits. 
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Appendix: Reflection Operators in Phase Space 



Among the several representations of quantum mechanics, the Weyl- Wigner represen- 
tation is the one that performs a decomposition of the operators that acts on the Hilbert 
space, on the basis formed by the set of unitary reflection operators. In this appendix we 
review the definition and some properties of this reflection operators. 

First of all we construct the family of unitary operators 

T q = exp(— ih~ l q.p), T p = exp(i^ _1 p.g), (57) 
and following 15|, we define the operator corresponding to a general translation in phase 



space by £ = (p, q) as 



T% = exp ( — £ A x 



TpT q exp 



-2h P - q 



exp 
= T q T p exp 



-(p.q-Q-P) 



2h m 



(58) 
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where naturally x = (p, q) . In other words, the order of T p and T q affects only the overall 
phase of the product, allowing us to define the translation as above. T£ is also known as a 
Heisenberg operator. Acting on the Hilbert space we have: 

Tt\q a >=efrfe» + % a + g> (59) 

and 

f € \ Pa >= e-i q{pa+ ^\p a +p> . (60) 

We, hence, verify their interpretation as translation operators in phase space. The group 
property is maintained within a phase factor: 

~ — % ~ — % 

T 6 T 6 = T 6+6 ex P^i A = T 6l+ b exp[— D 3 (^, f 2 )], (61) 

where D 3 is the symplectic area of the triangle determined by two of its sides. Evidently, 
the inverse of the unitary operator T^ 1 = = T_£ . 

The set of operators corresponding to phase space reflections R x about points x = (p, q) 
in phase space, is formally defined in [15J as the Fourier transform of the translation (or 
Heisenberg) operators 

R x = (Aixhy L J dt, e^ xA % (62) 
Their action on the coordinate and momentum bases are 

R x \q a ) = e 2 ^~^ h \2q - q a ) (63) 
R x \p a ) = e 2i{ *>-^' h \2p-p a ), (64) 

displaying the interpretation of these operators as reflections in phase space. Also, Using 
the coordinate representation of the coherent state fflBl and the action of reflection on the 
coordinate basis (16"3"|) . we can see that the action of the reflection operator R x on a coherent 
state \X) is the x reflected coherent state 

R x \X) = exp (^X A x^j \2x - X) . (65) 

This family of operators have the property that they are a decomposition of the unity 
(completeness relation) 

i = ^nf dxft " (66) 
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and also they are orthogonal in the sense that 



Tr 



R X1 Rx 2 = Z^h 8{x2 — xj 



(67) 



Hence, an operator A can be decomposed in terms of reflection operators as follows 

1 



.4 



2ixh 



dx A\y(x) R x . 



(68) 



With this decomposition, the operator A is mapped on a function Ay/(x) living in phase 
space, the so called Weyl-Wigner symbol of the operator. Using f|67|) it is easy to show that 
Ay/(x) can be obtained by performing the following trace operation 



A w (x) = Tr 



R x A 



Of course, as it is shown in [15(, the Weyl symbol also takes the usual expression in terms 
of matrix elements of A in coordinate representation 



A w (x) 



q 



Q 



A 



q+ 2 / exp \~h P ® ^ 



It was also shown in 
composition properties 



151 ] that reflection and translation operators have the following 



so that 



RxT^ 


= Rx-i/2^ 


1^ 

) 


T$R X 


— Rx+(,/2^ 


jrXAi 

1 


Rxi Rx2 


= T2( X2 -x 1 ) 











Now using (171j) and (1701) we can compose three reflections so that 

ODD _ lA 3 (x 2 ,Xl,x) S 
1 1 <X2 ± ^X 1 L Xl t> 1 h x 2 — X+Xl 



(69) 
(70) 
(71) 

(72) 
(73) 



where A 3 (x2,xi,x) = 2(x 2 — x) A (x\ — x) is the area of the oriented triangle whose sides 
are centered on the points x 2 ,Xi and x respectively (see figure 1). 
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FIG. 2: detail of the Scar Wigner function with TV = 223 constructed on the fixed point (0,0) for a non 
Bohr quantized value of 4>- 

(a) Exact result for the cat map. 

(b) We compare sections of the exact and the semiclassical Scar Wigner functions near the periodic point 
for q = 0.5 (horizontal section). The solid line represents the section of the exact Scar Wigner function, 
while in dashed lines we show the semiclassical approximation. 

(c) Idem (b) but for p = 0.5 (vertical section). 

(d) Semiclassical approximation Px<p( x )- 
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